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ly-j \ Abstract 

o 

We consider an impact of hadronic light-by-light scattering on the muonium 
\q '• hyperfine structure. A shift of the hyperfine interval Az/(Mu)hlbl is calculated 

with the light-by-light scattering approximated by exchange of pseudoscalar 
and pseudovector mesons. Constraints from the operator product expansion in 
QCD are used to fix parameters of the model similar to the one used earlier for 
the hadronic light-by-light scattering in calculations of the muon anomalous 
^ ■ magnetic moment. The pseudovector exchange is dominant in the resulting 

shift, Ai/(Mu)hlbl = -0.0065(10) Hz. Although the effect is tiny it is useful 
in understanding the level of hadronic uncertainties. 
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1 Introduction 



Pure leptonic objects, such as free electron and muon or leptonic bound systems, 
positronium and muonium, are of specific interest because they allow ah initio calcu- 
lations with a high accuracy. There is no effect of strong interactions in the leading 
terms and in a number of terms in perturbative series. Still hadronic effects enter 
through higher loops in electromagnetic and electroweak interactions. 

The most important leptonic property affected by hadronic effects is the anoma- 
lous magnetic moment of a muon, a M = (g^ — 2)/2, where the main hadronic con- 
tribution comes from the vacuum polarization (HVP), see Fig.[IJ The hadronic con- 




Figure 1: The hadronic vacuum polarization contribution to a M 

tribution is quite small Aa M (HVP) ~ 7 x 10~ 8 ~ 7 x 10 _5 a M , but nevertheless it is 
much larger than the experimental error in the measurement [2], 

a^ xp = 116 592 080(63) x 10~ n , (1) 

as well as the uncertainty of the QED calculations [1] and the electroweak contribu- 
tion. The HVP contribution is obtained with sufficient accuracy by applying data 
from e + e~ annihilation into hadrons. 

At this level of accuracy one need to take into account higher order hadronic 
effects and, in particular, the virtual light-by-light scattering (HLBL), see Fig.[2j In 




Figure 2: The light-by-light scattering contribution to 

contrast to HVP there is no direct experimental input for determining HLBL so one 
should rely on a theoretical model. 
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Two relevant theoretical parameters are the smallness of the chiral symmetry 
breaking, m\/m 2 p -C 1, and the large number of colors, N c ^> 1. The first parameter 
enters a powerlike, chiral enhancement for the charged pion loop in HLBL 

while the large N c limit implies dominance of meson exchanges, see Fig. [3], where 




Figure 3: Meson exchanges in the light-by-light scattering. Summation goes over inter- 
changes of photons and over C-even neutral mesons. 

mesons M include neutral pion and heavier C-even resonances. 

In a number of papers dwelt on the problem it was shown that the chirally 
enhanced two-pion contribution is significantly smaller than the color enhanced one 
[3-7]. The model for light-by-light scattering developed in [7] is based also on QCD 
constraints which follow from operator product expansion at large photon virtualities. 
Together with the neutral pion the exchange of pseudovector mesons plays major role 
in the model. 

In the present paper we consider an impact of the hadronic light-by-light scatter- 
ing on another 'pure leptonic' quantity, namely, to the muonium hyperfine splitting 
(HFS), see Fig.Hl The hadronic effects in muonium are of somewhat less practical 
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Figure 4: The hadronic light-by- light scattering in the electron- muon interaction 

importance since there has been no experimental progress for years [8]. However, 
the accuracy in the former experiment was limited by statistics due to low muon 
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flux. At present, better muon sources are available, e.g., at the Paul Scherrer Insti- 
tut, and more accurate results are in principle possible. To start preparation for a 
new experiment one has to clearly understand the ultimate limit of the theoretical 
accuracy. 

In principle, pure QED calculations are ab initio calculations and can be done 
with any accuracy (which does not mean that they can be done easy — see, reviews 
[9, 10] for the present status). However, the very involvement of the hadronic effects 
sets a certain limit of accuracy. As well as in the case of a M one has to calculate the 
HVP contribution in the leading order [11,12], see Fig.[5j and the next-to-leading 



Figure 5: The characteristic diagram for the hadronic vacuum polarization contribution 
to the muonium HFS interval 

term [12,13]. 

The HVP contribution can be found from experimental data on the e + e~ anni- 
hilation into hadrons. The HLBL contribution is of the same order as the next-to- 
leading HVP contributions [12] and cannot be derived from existing scattering and 
annihilation data. So we extend the model of Ref. [7] for the hadronic light-by- light 
scattering to apply it to the muonium HFS. 

An interesting feature of this application is that the dominant contribution comes 
from the "vertical" exchange by pseudovector mesons, see Fig. [61 The reason for this 



Figure 6: The "vertical" exchange of pseudovector meson a in e/i scattering 
dominance is that the pseudovector exchange shown in Fig. [6] is the most relevant one 
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for the spin-spin interaction of the electron and muon which determines the HFS. By 
contrast, a similar exchange of a neutral pion vanishes in the scattering amplitude 
when the electron and muon are at rest. 

The pion and pseudovector cross-channel "horizontal" exchange, see Fig. [7J are 

U XT' 

1_ cL l pi 

Figure 7: The "horizontal" exchange of pseudoscalar tt and pseudovector mesons a in the 
efj, scattering 

also accounted in the model. This contribution is numerically smaller than the 
"vertical" one. Thus, the situation in HFS is opposite to that for the a M , where the 
pseudoscalar exchanges dominate. 

Another interesting point is that the chirally enhanced charge pion loop in the 
blob of Fig. H] does not contribute to spin-dependent part of the scattering amplitude. 
Indeed, the quantum numbers J p of exchange should be 1 + as for pseudovector 
mesons. However, such quantum numbers are not allowed for the pair 7r + 7r~. Thus, 
in contrast to a M , an ambiguous charged pion loop does not enter the muonium HFS. 

In the next section we introduce general expressions for the HLBL effect in HFS. 
In Sec. 3 we present calculations of the pseudovector exchange, and in Sec. 4 we 
consider the pseudoscalar exchange. In the last section we summarize the results. 

2 Generalities 

Let us start with some general expressions. The muonium HFS is determined by the 
spin-dependent part of the forward e~fi + scattering in the low-velocity limit. It is 
convenient to start with the e~fi~ amplitude and then make the charge conjugation 
for the muon. 

The spin-dependent part of the forward e~ (p) + \i~ (r) — > e~ (p) + fi~ (r) scattering 
associated with HLBL can be presented as 

M spin (e- \i~ — > e~ pT) — A u (e) 7 ff 75 u (e) u M ^ 5 u M — ► -Am e m^A a e , (2) 

where u^ e \ vS™ are Dirac spinors describing the electron and muon (we are using 
relativistic normalization and units h = c = 1) and we took the nonrelativistic limit 
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in the last expression. The transition from yT to /i + does not change the result 
because of the positive C-parity of the axial current u^j^jbU^. 

The above amplitude leads to the following addition in the e~n + Hamiltonian, 



AH UFS = A6 3 (r)a e a,. (3) 
This should be compared with the leading term for the s-wave HFS Hamiltonian, 

H» FS = ^^6 3 (r)a e a„ (4) 
which gives for the HFS interval (Fermi energy Ep) 

hu HFS = Ep = ( m » V ~ h ■ 4.459 x 10 9 Hz . (5) 

3 m M ynie + m^J 

The shift in the splitting due to HLBL is 

/^hlbl = A£ H lbl = E F A . (6) 

lira 

The amplitude A is defined by the diagram in Fig.HJ 

_ Ao^ f d A kd A q M^, u ,(k,q) , s , p , vS , / _ rjr^ \ 

3 J (27T)6 (P)3g2[( g 2 ) 2_ 4(rg) 2] e K P £ <V \?8S> ^ ) • UJ 

Here M^^^^k, q) is the the amplitude of the forward scattering of two virtual pho- 
tons, 

t*(a fc ) + t*(^ q) -> 7*(/A fe) + 7*(V, g) , (8) 

(we mark their polarization indices and momenta), r = {m M , 0} is the 4-momentum 
of the muon at rest, and we neglected by the electron mass. In many cases, in 
particular for pseudovector exchanges, one can neglect by the muon mass as well. 
Then the expression in Eq. ([7]) simplifies further, 



A = a' 



d 4 kd 4 q M 



^e^V^W (9) 



(2tt) 6 (k 2 f{q 2 f 

3 Pseudovector exchange 

To calculate the "vertical" pseudovector exchange, see Fig.[6j let us start by intro- 
ducing the effective vertex for lepton interaction with the pseudovector meson a, 

hadplYlsl, I = e, (J,, (10) 



5 



where h a is the coupling constant and a p is the polarization of the axial meson. 
Implying the same h a for the electron and muon (corrections due to their mass 
difference are small and can be accounted for ) we get the pseudovector contribution 
to the amplitude A of the forward e/i scattering, see Eq. (T2]), 

h 2 

4 pv _ nn 

where m a is the pseudovector meson mass. Note the negative sign which follows 
from unitarity (see Ref. [14] for a detailed discussion of the sign). Hence even before 
explicit calculation we know that the pseudovector exchange correction to HFS is 
negative. 



3.1 Coupling of pseudovector mesons to photons and leptons 

The next step is to calculate h a . To fix the 07*7* vertex one can use that at large 
virtualities the operator product expansion relates product of two electromagnetic 
currents to the axial current [15], see also [7], 

|dWye-^^T{j^^ (12) 

Here 

f 5 =qQ 2 Y l5 q (13) 

is the axial current, where different flavors enter with weights proportional to squares 
of their electric charges, <?3 = gi + q2 and q = (qi — 92) /2 ~ qi ~ — <?2 • 

The 07*7* vertex which satisfies this constraint at large q 2 and regular at small 
q can be chosen in the form 



v CiU.V 



{q\ -ml){ql-rnl) 



qh^s P q[ + (qi <-»• 92, n <-+ v) ■ (14) 



Here (a|jg (0)|0) is the matrix element between vacuum and the outgoing axial meson 
with 4-momentum g 3 = q\ + q 2 and polarization a p (photon momenta q\ and q% 
momenta are taken as incoming). The form factor parameter m v is the mass of the 
appropriate vector meson. Of course, this form of the vertex is model-dependent. 
This refers not only to the above form of q\,q\ dependence but also to choosing a 
particular structure, one of three possible structures for the vertex. The choice ffT4|) 
picks up the structure which survives in asymptotics. 
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The lepton interaction with a can be calculated then from the triangle diagram 
(the upper and lower blocks in Fig.[6]). Taken all external momenta to be vanishing 
and neglecting by lepton mass we get 

V = -2eMMW)^^^^h»^l 

= -^<«K(0)|0)'7p75(- (15) 

It gives the result for the pseudovector coupling to leptons in terms of the vector mass 
m v and the matrix element of the axial current between the vacuum and pseudovector 
meson, 

Actually there are three electrically neutral pseudovector mesons which differ 
by their features under flavor SU(3). Therefore it is convenient to present the axial 
current j$ p = <?<5 2 7 P 75 q as a linear combination of axial currents with the same SU(3) 
quantum numbers as the mesons a*''. In particular, we can introduce the isovector, 
dtp = qhlplsq, hypercharge, jj$ = qXslplsq, and the SU(3) singlet, j$ = c?7 P 7 5 g, 
and write 

E Tr[A fc y ] .(ft) / 1R s 

fc=3,8,0 ir ^^ 

where Ao is the unity matrix. Accounting for mixing of the hypercharge and singlet 
pseudovector mesons is simply done by substituting Ag and Ao by appropriate linear 
combinations. 

Thus, we get for the coupling ha of the meson to leptons (see Eq. ( jTUi) for 
definition) 

{k) _ Tr[A fc Q 2 ] 

ml Tr[A^] L ' (17j 

where fa is the coupling of the meson to the corresponding axial current, 

(« (fc) |ji?|0) = /i fc) ap. (18) 



3.2 Coupling of pseudovector mesons to axial currents 

The value of fa can be fixed from consideration of the transition of the axial current, 

(k) 

j$ p , into two photons. We consider a special kinematics when one of those photons is 
soft with momentum k — > and polarization e u and another is virtual, carrying the 
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same momentum q as the axial current. The transition amplitude can be represented 

as 



rp(k) I 



i(0\ld"ze^T{fJ(z)ej ll (0)}\ 1 ). 



(19) 



Generically, as it is shown in [16], the transition Tp^l can be written in terms of two 
Lorentz invariant amplitudes, w^\q 2 ) and w!^\q 2 ), 



rp(k) V 

pp,v 



where /. 



ftp 



ie 2 N c Ti{\ k Q 2 
4vr 2 



w ( L k \q 2 ) q P qJa, + w ( T K \q z ) [-q 2 U P +q»qJa P -q p qJ ( 



,( fc V„2 



/ tTfl 



(20) 



In perturbation theory, w^ T (q 2 ) are computed from triangle diagrams with two 
vector currents and an axial current. For massless quarks, we have 



w[ k \q 2 ) = 2<V) = ~ ■ 



(21) 



An appearance of the longitudinal part for the axial current which classically is 
conserved is a signal of the famous Adler-Bell-Jackiw axial anomaly [17]. The pole 
at q 2 = in w L \q 2 ) is associated with propagation of massless Goldstone particles, 
the pion in case of w^. 

There is no perturbative corrections to these functions in the chiral limit. More- 
over, the longitudinal functions w L ' protected even against nonperturbative correc- 

(k) 

tions. It is not the case for transversal functions Wj, where the pole should be shifted 
from zero to vector and pseudovector masses. A particular model which account for 
this shift suggested in [16] has the form (in the chiral limit), 



Wj. (q 



1 



mt 



ml — q 2 



(22) 



where m a>v denote masses of pseudovector and vector mesons in the given channel k. 
Equation ffl9l implies the following expression for the residue of the pole at 



lim (g 2 - m 2 a ) T« <?= f^V mv { qi = q,q 2 = k)e» . 



(23) 



Comparing this with the residue from Eqs.( l20i) and ( 1221) we get the result for 
In particular for we have 



J a 



(3) 



N c m 4 p 
2tt 2 



(24) 
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(3) 

An independent way to find /„ is to use Weinberg's sum rules to relate it with the 
p coupling to electromagnetic current, (plj/JO) = (m 2 /g p )p M , 



J a 



(3) 



-i 2 



2m 2 
9 P 



(25) 



Then Eq. (1241) implies an interesting relation 



El 

Atx 



2tt 
N~ r 



(26) 



reasonably good phenomenologically This can be also compared with the QCD 
sum rule result [18], g 2 /(47r) = 27r/e, where e, the base of natural logarithm, enters 
instead of N c - a pretty good approximation for N c = 3. 



3.3 Pseudovector exchange results 

Combining Eqs. ( [IT]) . ( fTTj) and ( 1251) we get for the exchange by the isovector ai(1260) 
meson, 

3 a 4 1 



A ai 

vert 



7T 2 W? 



(27) 



ai 



For the isoscalar pseudovector mesons /i(1285) and / 1 *(1420) we assume the ideal 
mixing, similar to 10 and <fi. It means that f\ has the (uu + dd) / a/2 structure and 
/* is ss; this assumption is consistent with experimental data for decays of these 
resonances. In terms of the relevant axial currents the linear combinations (Ao/3) — 
(As/\/3) and (2A /3) + (X 8 /^/3) enter correspondingly. Then, similarly to the a\ 
exchange, we get 



vert 



25 a 4 



24 7r 2 m 2 h 



A 



n 

vert 



1 a 4 1 



12 7T 2 mi 
Ji 



Altogether the "vertical" exchange by pseudovector mesons produces 



(28) 



PV 



71 



m P m, 



Ol 

a \ 3 m e m. 



71 



Ol 



9 25 mg 1 mg t 
16 16 m 2 * 8 m 2 ,* 

E F ■ 2.16 = h ■ (-0.0041 Hz 



(29) 



Now let us add up the "horizontal" pseudovector exchanges shown in Fig. [71 In the 
limit of heavy pseudovector mass this exchange differs from the "vertical" one just by 
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averaging over angles and constitutes 1/3 of the "vertical" exchange. Accounting for 
the finite pseudovector mass we found an extra to 1/3 suppression of the "horizontal" 
exchange by the factor 0.614, 

pv /a\ 3 m P m u 0.614 , , . 

A£™ - f \ _^_J1 Ef .2.16 — = h- -0.000 84 Hz . 30 

\7T/ mi, 3 

Thus, the total for pseudovector exchange is 

A £ PV = - ( «V !^ Ep -2.Q = h- (-0.0049 Hz) . (31) 
\7T/ m„, 

Note that we limit ourselves by exchanges of pseudovectors with the lowest mass 
in each flavor channel. Exchanges by higher 1 + excitation contribute in the same 
direction but probably are numerically suppressed. 



4 Pseudoscalar exchange 

Let us start with the 7r°7*7* vertex, 

V[lll Cftyy F^ry* ry* ( k ,Q ) ^ fJ,U pa ^ Q • (^2) 

Here k and q are photon momenta, fi and v are their polarization indices, the constant 
c 7r77 is fixed by the width of 7r° — > 77 decay and F^*^*(k 2 , q 2 ) is the form factor of 
the transition, F 7r7 » 7 »(0,0) = 1. Theoretical expression for 

follows from the Adler-Bell-Jackiw anomaly [17]. Indeed, it could be read off from 
the residue of the pole in the longitudinal part in Eq. (1201) . 

The pion exchange gives then the following expression for the forward scattering 
of two virtual photons: 

M-V = c^[K rr (k 2 , q 2 )f ^fy { ^/ + (a* - A*' ■ (34) 

The third permutation involving k <-> — q vanishes for the forward scattering. It 
means an absence of the "vertical" exchange for the pion mentioned earlier once 
atomic momenta are neg lectedQ 



1 The higher order contributions to muonium HFS due to the "vertical" pion exchange together 
with additional photon are suppressed by extra small factors such as a and m e /m n . 
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Now we have to substitute M v p ° p , v , to Eq. (jTJ) and integrate over k and q. By 
power counting at large momenta it is simple to see that in absence of the form 
factor F 7r7 . 7 *(/c 2 , q 2 ) the integral logarithmically diverges. The form factor provides a 
convergence above momenta of order of m p , while its infrared convergence is regulated 
by pion and muon masses. The \n(m p /m n ) term can be determined analytically, 



For numerical estimates we use = 135 MeV, m p = 775 MeV, F n = 92 MeV. 
The logarithm is not that big, \n(mp/m n ) = 1.75 so a numerical integration with a 
certain model for the form factor is needed. 
For the form factor 

4 

771 

i^ 7r7 * 7 *(fc q ) = - - — — - — (36) 
[m l — k 2 )(m 2 p — q 2 ) 

numerical integration gives 



The suppression of the logarithmic result can be approximated by substitution 



in Eq. (jggjl. 

The result (1371) can be compared with the earlier calculation of the pion contri- 
bution by Faustov and Martynenko [19]. They found AE n = h ■ (+0.0011 Hz) for 
the same form factor (1361) . While the magnitude is close we differ in the sign. Note 
that Faustov and Martynenko also considered change of HFS due to effect of HLBL 
pion exchange on a p , the effect we are not considering. 

Strictly speaking the form factor (1361) violates the QCD constraints. It follows 
from the OPE expansion (I12p that at q 2 = k 2 the form factor should decrease as 
1/q 2 at large Euclidean momenta, not as l/q 4 as in Eq. (136]) . So we made numerical 
integration with the form factor which satisfies the above mentioned constraint as 
well as other theoretical and experimental limitations [5] (see also [7]), 




(35) 




(37) 




(38) 



{k , q ) 



m 



, p M* - (4tt 2 F 2 /N c ) q 2 q 2 (q 2 + q 2 ) + h 2 ■ q\q 2 + h 5 • (q 2 + q 2 ) 



{q\ - m 2 p )(q 2 - M 2 )(q 2 - m 2 p )(q 2 - M 2 ) 



(39) 
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where M 2 = 1465 MeV , h 5 = 6.93 GeV 4 , h 2 = -10 GeV 2 . The result of integration 
turns out to be very close to the one in Eq. (1371) . the difference is insignificant. 

Calculations for the other pseudoscalars, 7/(547) and 7/(958), can be done in a 
similar fashion. We use their experimental two-photon width to determine the two- 
photon couplings and simple vector-dominance form factors for off-shell photons with 
m p = muJ = 775 MeV and = 1020 MeV, 



F vrr (k 2 ,q 2 ) 



5 rrf 2 77^ 

3 (m 2 — k 2 )(m 2 — q 2 ) 3 (m 2 — /c 2 )(m 2 — q 2 ) ' 



F (k 2 a 2 ) - m " , I m l (m 

{ ' q ' 6 (m 2 - k 2 ) (m 2 - q 2 ) 6 (m 2 - k 2 )(m 2 - q 2 ) ' { ' 

based on the octet and singlet quark structure of rj and rj'. The results of numerical 
integration are 

AE„ = -[-) , e ^ Ep ■ 0.063 = h ■ (-0.000 14 Hz) , 

\7r/ {A-nF^y 

/ qj \ 3 ^2 /7T, 

A£L, = - - —±-JL- E F - 0.046 = h ■ (-0.000 10 Hz) . (41) 
\7r/ (47riv) 2 

This can be compared with calculations by Faustov and Martynenko [19], they ob- 
tained 0.0002 Hz for rj and 0.0001 Hz for rj'. Again, we have a sign difference. 
Thus, the total for pseudoscalar exchanges is 

AE PS = - (-) 3 r ? e ^ 2 E F • 0.72 = h ■ (-0.0016 Hz) . (42) 



5 Summary 



Collecting the results (1311) . (1421) for pseudovector and pseudoscalar exchanges we 
come to 



A^hlbl = - (^-J rn e m^ E F 



2.6 0.72 
+ 



_m 2 ai (4ttF^) 

= /i- (-0.0049 Hz -0.0016 Hz) = h ■ (-0.0065 Hz). (43) 

The main contribution is due to pseudovector exchange, the "vertical" one in Fig. [61 
It is fivefold larger than "horizontal" pseudovector exchange and threefold larger 
than the "horizontal" pion exchange. 
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What is the accuracy of the result? We mentioned in Introduction an absence 
of the charged pion loop associated with chiral enhancement. This makes the result 
more reliable. Looking on variations of parameters such as coupling of the pseu- 
dovectors to axial currents we would estimate the uncertainty of the model for the 
dominant pseudovector exchange as 10%. Staying on the conservative side we ascribe 
a total uncertainty to be about 25% of the pseudovector "vertical" contribution, i.e., 
about 0.001 Hz. Thus, our final result is 

/iAz/hlbl = AEhlbl = -0.0065(10) Hz . 

We see that the HLBL correction is tiny and rather unobservable. However, it 
shows the level of limitations on theoretical accuracy which comes from hadrons. In 
our study we also obtained a few relations for couplings of pseudovector mesons to 
photons, leptons and axial currents which can be applied to variety of processes. 
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